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Ell  ABSTRACT 

Thomas  Young  derived  in  1807  a  formula  which  associated  the  applied  load,  the 
Euler  critical  load,  the  initial  bow,  and  the  elastic  deformation  under  load  for  a 
column.  Westergaard  in  1922  used  the  Lagrangian  law  to  derive  Young's  expression  in  a 
much  wider  context,  and  in  so  doing  demonstrated  its  general  validity  in  all  linear 
instability  problems.  This  relationship  has  great  significance  in  the  correlation  of 
experimental  daui  with  linear  theory.  This  was  realized  first  by  Ayrton  and  Perry  in 
1887.  In  1932  Southwell  Independently  redeveloped  the  technique  and  used  it  to 
discuss  the  classic  column  tests  of  von  Karman. 

This  report  coordinates  and  extends  the  research  made  in  this  field  during  the 
last  century  and  a  ha]f .  It  demonstrates  the  power  and  versatility  of  the  method 
of  expressing  the  given  hyp  .•"bolic  equation  in  a  linear  form  -  by  change  of 
variables  -  and  determine.)  the  classic  critical  load  as  the  slope  of  this  line.  The 
results  show  that  there  can  be  no  doubt  as  to  the  applicability  of  the  method  to  all 
linear  instability  proolems  of  columns  and  plates.  Thus,  the  report  challenges 
the  statements  of  previous  researchers  who  have  been  hampered  in  their  work  by  the 
lack  of  valid  expert  nents,  arising  from  the  coarseness  of  measurement  techniques. 
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SUMMARY 


Thomas  Young,  medical  doctor,  philosopher,  and  Professor  of  Mechanical  Arts 
at  the  University  of  Cambridge,  derived  in  1807  a  very  simple  but  powerful 
formula  which  associated  the  applied  load,  the  Euler  critical  lead,  the 
initial  bow,  and  the  elastic  deformation  under  load  for  a  column.  His 
hyperbolic  law  is  in  essence  we  =  (wqPe/(P  -  PE)).  Westergaard  in  1922 
used  the  Lagrangian  law  to  derive  this  expression  in  a  much  wider  context, 
and  in  so  doing  demonstrated  its  general  validity  in  all  linear  instability 
problems . 

This  relationship  has  great  significance  in  the  correlation  of  experimental 
data  with  linear  theory.  This  was  realized  first  by  Ayrton  and  Perry,  who 
in  1887  published  a  paper  on  struts  in  which  they  gave  a  graphical  inter¬ 
pretation  of  the  fundamental  equation,  and  by  so  doing  were  able  to  verify 
that  the  Euler  load  was  correct  for  a  perfect  column. 

Their  important  work  was  virtually  ignored.  However,  in  1932  Southwell 
independently  redeveloped  the  technique  and  used  it  to  discuss  the  classic 
column  tests  of  von  Karman.  In  his  report,  Southwell  stated  without  proof 
the  axiom  which  was  so  clear  in  Westergaard 1 s  analysis;  viz,  that  the 
process  should  have  wider  use. 

This  report  reviews  these  several  contributions  to  our  knowledge  together 
with  the  numerous  papers  which  have  followed  them.  A  wide  range  of  problems 
is  illustrated  both  experimentally  and  analytically.  In  those  cases  for 
which  the  appropriate  theory  or  the  valid  experiment  exists,  the  current 
results  have  been  presented;  but  where  gaps  in  reported  information  exist, 
the  necessary  analysis  or  tests  were  conducted  and  presented  herein. 

It  is  significant  to  record  that  the  work  covers  not  only  normal  load 
environments  but  also  instabilities  under  thermal  conditions. 

Thus,  we  may  state  that  the  report  coordinates  and  extends  the  research 
of  some  twenty  engineers  made  in  a  century  and  a  half.  It  demonstrates 
the  power  and  versatility  of  the  method  of  expressing  the  given  hyper¬ 
bolic  equation  in  a  linear  form  -  by  change  of  variables  -  and  of  deter¬ 
mining  the  classic  critical  load  as  the  slope  of  this  line.  The  results 
obtained  show  that  there  can  be  no  doubt  as  to  the  applicability  of  the 
method  to  all  linear  instability  problems  of  plates.  Thus,  the  report 
corrects  the  statement  of  previous  researchers  who  have  been  hampered  in 
their  work  by  the  lack  of  valid  experiments,  arising  from  the  coarseness 
of  measurement  techniques.  There  is  clear  indication  that  continued 
study  will  open  the  door  to  a  clearer  understanding  of  the  perplexing 
questions  which  arise  in  compressive  elements  in  which  curvature  is 
present. 
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The  subject  of  elastic  stability  is  more  than  200  years  old.  From  the 
scientific  point  of  view  it  begsa,  in  17^,when  Leonard  Euler-*-  made  the 
now  classic  analysis  of  the  behavior  of  an  initially  straight,  perfectly 
elastic,  centrally  compassed  strut.  Despite  his  analytical  achievement, 
there  is  little  doubt  that  in  the  century  which  followed,  practicing 
engineers  still  had  to  rely  upon  experience  and  enpiricism  bora  in  experi¬ 
mentation. 

There  is  no  doubt,  however,  that  researchers  who  were  active  in  the 

field  studied  Euler’s  work  very  thoroughly.  Indeed,  Eodgkinson,  a  foremost 

experimentalist,  remarked  in  1842  with  reference  to  Euler.*"5 

"It  appeared  from  the  researches  of  this  great  analyst,  that  a  pillar 
of  any  given  dimensions  and  description  of  material  required  a  certain 
weight  to  bend  it,  even  in  the  slightest  degree;  and  with  less  than 
this  weight  it  would  not  be  bent  at  all.  Lagrange,  in  an  elaborate 
essay  in  the  same  work,  arrives  at  the  same  conclusion.  The  theoiy 
as  deduced  from  this  conclusion  is  very  beautiful,  and  Poisson's 
exposition  of  it,  in  his  Mecanique,  2nd  edition,  vol  i,  will  well 
repay  the  labour  of  a  perusal. 

I  have  many  times  sought,  experimentally,  with  great  care  for  the 
weight,  producing  incipient  flexure,  according  to  the  theory  of  Euler; 
but  have  hitherto  been  unsuccessful.  So  far  as  I  can  see,  flexure 
commences  with  weights  far  below  those  with  which  pillars  are  usually 
loaded  in  practice.  It  seems  to  be  produced  by  weights  much  smaller 
than  are  sufficient  to  render  it  capable  of  being  measured." 

At  the  time  that  Hodgkinson  made  this  remark  a  large  quantity  of  experimental 
data  was  being  obtained  by  many  workers.  These  r -perimental  data  were, 
however,  at  variance  with  the  predictions  of  Euler,  and  analysts  thus  came 
to  regard  his  theory  as  incomplete.  Now  analysis  according  to  the  classic 
theory  indicates  that  struts  of  the  type  already  described  will  remain 
straight  until  a  certain  critical  load  level  is  reached;  then  the  stnrts 
may  suddenly  bend  into  the  form  of  a  single  bow,  a  shape  which  will  o- 
maintained  with  increasing  load  although  the  amplitude  of  the  deformation 
will  increase. 

In  reality,  however,  struts  tested  under  axial  compression  have  lateral 
deformation  from  the  very  onset  of  loading.  It  is  found  in  practice  that 
the  load  deformation  curves  for  such  bodies  have  the  characteristic  shapes 
given  in  Figure  1.  When  they  are  long,  nhe  curve  tends  to  become  asymptotic 
to  the  Euler  curve,  but  when  they  a.re  shorn,  it  does  non  do  so.  The  curve 
for  short  struts  is  characterized  by  maximum  values  of  load,  a*"  is  clearly 
seen  in  the  figure.  In  is  not  difficult,  therefore,  to  arrive  at  the  point 
of  view  that  agreement  between  the  critical  load  determined  by  Euler  consid¬ 
erations  and  the  buckling  load  as  determined  by  experiment  is  a  chance 
occurrence  to  be  expected  only  in  the  case  of  very  long  struts.  Thus,  one 
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is  forced  to  ask  tte  question:  When  experimental  observations  give  curves 
of  the  tjpe  presented  in  Figure  1  and  the  critical  load  is  never  attained, 
can  any  valid  comparison  of  theory  and  experiment  be  made? 

The  first  researchers  to  provide  a  substantial  answer  to  this  question  were 
Ayrton  and  Perry.-*  In  1886  these  two  investigators  published  an  instructive 
paper  on  struts.  It  is  interesting  tc  note  that  when  they  performed  this 
work,  the  situation  with  regard  to  knowledge  on  this  fundamental  question 
was  such  that  they  comnentei:  "We  think  we  are  right  in  saying  that  no 
subject  connected  with  the  strength  of  structures  is  at  present  regarded 
as  being  in  a  more  vague  and  unscientific  state  than  the  subject  of  which 
we  have  undertaken  the  study". 

As  a  result  of  their  study,  Ayrton  and  Perry'*  were  among  the  first  to 
realize  that  for  a  complete  theory  it  was  essential  to  introduce  imper¬ 
fections  of  form  and  loading.  Their  own  words  express  this  particular 
point  and  their  manner  of  dealing  with  it:  "The  conclusion  at  which  we 
have  arrived  is,  then,  that  any  want  of  axiality  in  placing  the  load,  or 
want  of  straightness  in  the  unloaded  strut,  or  want  of  homogeneity  in  the 
material  may  be  allowed  for  by  a  terra  C  such  that  it  may  be  taken  as  the 
initial  deflection  of  a  homogeneous,  carefully  loaded  strut,  and  then, 

y,  =  s  (i) 

i  -  tegp 

EUT2 

where 

y  =  the  ordinate  of  the  strut  center  line  in  the  displaced 
configuration 

E  =  Young's  modulus 

I  =  second  moment  of  cross-sectional  area 
L  =  pin  center  length 
P  =  axial  load." 

From  this  expression,  they  deduced  that  the  measured  deflection  of  the 
midpoint  from  its  initial  value  was  given  by 


where  P  was  the  end  load  and  P,  was  Euler  load. 
Equation  (2)  can  be  written  as 


If  we  choose  new  variables, 

1  .  1 
x  =  g  and  y  =  ^ 

then  we  see  that  equation  (3)  is  a  linear  equation  in  these  variables.  Thus 
a  plot  of 

11  C 

-g  and  |  is  a  straight  line  whose  slope  is  p  . 

This  analytical  deduction  was  verified  from  tests  on  Dantzic  0a*;  and  wrought 
iron.  The  appropriate  data  are  given  in  Figures  2  to  5 •  The  plots  of 


are  seen  to  be  straight  lines.  Ayrton  and  Perry's  result  of  equation  (2) 
can  also  be  expressed  as 

P,  ^  -  6  =  C  =  effective  initial  deflection  (4) 


Thus,  we  may  choose  a  new  variable  x,  such  that 


Hence,  equation  (4)  may  be  written  as 

P,x  -  6  =  C  =  constant  (5) 

This  is  a  linear  equation  in  x  and  6.  The  inverse  of  the  slope  of  the  line 
is  clearly  the  Euler  critical  load,  and  the  intercept  w  ith  the  6  axis  is  the 
effective  initial  deflection.  This  is  shown  in  Figure:  6  and  7 . 

It  is  strange  that  this  equation,  which  gives  a  very  elegant  method  of 
correlating  theory  and  practice,  was  virtually  ignored.  Indeed,  almost 
half  a  century  passed  before  the  method  was  used  to  any  extent.  Then  it 
was  independently  rediscovered  and  presented  in  the  second  form  by  R.  V. 
Southwell, °  who  published  his  work  in  the  Proceedings  of  the  Royal  Society. 
Today  the  method  of  presentation  is  always  referred  to  as  a  Southwell  plot. 
Since  it  is  this  latter  work  to  which  all  reference  is  made,  the  theoretical 
justification  for  the  formula  given  in  equation  (2)  will  be  derived  in  the 
Southwell  manner. 

Consider  the  strut  shown  in  Figure  8 •  Let  y  denote  the  initial  deflection, 
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Plot  Of  l/p 


Figure  Stress  Deflection  Curve.  Theoretical  Calculations  of 
Ayrton  and  Perry 5  on  Wrought  Iron  Strut. 


Plot  of  l/p  vs  l/6 


and  let  y  denote  the  final  deflection  of  the  central  line  from  the  line  of 
thrust.  Let  P  he  the  axial  load  and  let  El  he  the  flexural  rigidity  of  the 
body.  Then  the  condition  of  equilibrium  of  the  bent  body  is 

El(y"  -  f'o)  +  Py  =  0  (6) 

which  may  be  written  as 

y"  +  c^y  =  y"o  (7) 

Now  the  quantities  y  and  y  are  both  regarded  as  functions  of  x.  Thus, 
the  form  of  one  is  dependent  on  the  other.  If  y  vanishes  at  either  end 
of  the  strut,  then  a  gene  rail  solution  of  equation  (7)  can  be  obtained  if 
y  and  y  are  represented  by  the  following  Fourier  series: 


[0nsinSr] 


n=l 


r  =  y  [c 

o  L  L  n 


'(?  sinJSHl 


L  J 


n=l 


We  find,  then,  that 


(8) 


(9) 


CnV^ 

n  tt 
"th 

Thus,  if  P^  is  the  n  critical  load  given  by  the  Euler  conditions, 


(10) 


Jn  1 


C  .  P 
n  1"P 


'  n 


(11) 


Now  the  increase  in  deflection  at  the  center  of  the  strut  (elastic 
deflection)  is 


d  =  +  etc. 


(12) 


and  if  P  is  a  fairly  considerable  portion  of  P^,  this  may  be  written  as 
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d  -  1}  =  Constant  (lk) 

This  equation  is,  of  course,  identical  with  that  of  Ayrton  and  Perry  and 
is  in  fact  that  of  a  rectangular  hyperbola.  However,  whereas  Ayrton  and 
Perry  displayed  it  in  a  linear  form  by  plotting  the  reciprocals  of  the  load 
and  the  displacement,  Southwell  chose  to  plot  the  ratio  of  the  displacement 
to  the  load  against  the  displacement.  There  is  a  third  way  in  which  this 
equation  can  be  expressed  in  a  linear  manner.  In  this  alternative,  the 
load  is  plotted  against  the  ratio  of  the  load  to  the  displacement.  This 
latter  form  was  given  by  Donnell  in  his  paper  on  the  Southwell  method. 7 

The  differences  among  the  Ayrton  and  Perry,  the  Southwell,  and  the  Donnell 
representations  are  clearly  shown  in  Figure  9. 

Just  as  his  predecessors  had  chosen  the  most  accepted  test  data  of  their 
day  to  try  out  their  approach,  so  did  Southwell  pick  the  well  acclaimed 
data  of  von  Karman°  (Table  I) . 

The  results  are  shown  in  Figure  10.  In  all  cases  the  straight  lines  are 
good.  An  analysis  of  the  critical  load,  as  derived  from  the  plot,  and 
of  the  value  given  by  the  theoretical  formula  3s  given  in  Table  II.  The 
correlation  is  excellent. 


One  condition  prescribed  in  the  derivation  of  the  fonr'J 0  was  that  the  y 
should  be  zero  at  both  ends  of  the  strut.  Thus,  since  this  case  is  likely 
to  be  rare,  it  is  of  clear  interest  to  know  whether  or  not  in  practice 
the  method  would  still  work  for  eccentric  loadings.  R.  V.  Southwexl° 
studied  the  problem  and  confirmed  that  small  eccentricities  of  loading 
point  could  be  accommodated.  The  experimental  data  which  he  used  for 
this  purpose  were  obtained  by  Robertson. 

6  9 

After  the  publication  of  the  Southwell  paper,  Lundquist  reexamined  the 
question  and  generalized  the  result. 


The  formula  which  he  developed  is  as  follows: 


a^sin 


TTX 

L 


P  -  P 
crit  3. 


(15) 
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6  =  elastic  deflection 
P  =  applied  load 
P-j=  Euler  load 
c  -  constant 


i^LTth^i*0  i  fpproxlmation/  t0  the  load 

it  and  the  Euler  loaa  (P^  by  the  expression 


(P)  which  produces 


6  = 


Three  possible  arrangements  of  this  form 
Ayrton  and  Perry 

p 

formula  1/6  =  _  ij  i 

variables  1/6  and  l/p 

are: 

l 

P 

! 

y 

fU 

Ax  „  p 

Ax  _  __1 
Ay  =  c 

K-  c  -J 

- - -  Wi 

Southwell 

fP-,  1 

formula  6  j^-— ;  _  i j  _  c 

variables  6/p  and  6 

6 

P 

i 

"n 

H 

a* 

ii 

s 

1-  c  -1 

Donnell 

f ormula  P  -  r  =  c  ~ 

j.  o 

variables  P  and  P/6 


Figuze  9.  Linear  Graphical  Relationships  for  Columns. 
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XAfcMjii  i.  T.  vow  KARMM'S  STRUTS  (Ref.  o),  Nos  vl,  2 ,  3a,  3b,  4a,  4b,  5,  and 
6.  (MELD  STEEL:  E  =  2,  170.000  kg/cm2) 


P  =  End  Load  in 
Kilograms 


X  =  Measured  Deflection 
in  Millimeters 


V  =  x/P  x  10£ 


Strut  No.  1— 
2260 
3020 
3170 
3320 
3470 
3620 

Strut  No.  2 - 

4520 

4830 

5130 

5280 

5430 

Strut  No.  3a — 
6030 
7540 
8290 

8520 _ 

Strut  No.  3b— 
*7840 
81 40 

8290 

8445 

8600 

Strut  No.  4a--- 
*9050 
*9660 
10260 
10560 
10710 
10860 
11010 
11160 

Strut  No.  4b — 
*3020 
*4530 
*6030 
*7540 
*8300 
9050 

9805 

9960 

10110 

10260 

10410 


0.01 

0.025 

0.04 

0.06 

0.09 

0.25 


0.02 

0.025 

0.03 

0.07 

0.10 

0.13 

0.25 

0.73 


4.43 

8.28 

12.62 

13.07 

25.94 

69.06 

4.43 
10.35 
21.44 

45.46 

158.38 

1.66 

3-98 

13.27 

61.03 

2.55 

6.14 

8.44 

13.03 

24.42 

2.21 

2.59 

2.92 

6.63 

9o4 

11.97 

22.71 

65.41 

9.93 

11.04 

11.51 

11.94 

14 .46 
16.58 
23. 4o 

70.10 
.’’0.6  i 
32 .16 
39.  J9 


TABLE  I  -  Continued 


10560 

10710 

10860 

- 

0.71 

1.46 

49.24 

66.29 

134.44 

[Strut  No.  5—  ! 

*9050 

0.01 

1.105 

*10560 

0.03 

2.84 

10860 

0.05 

4.67 

IU60 

0.07 

6.27 

11^70 

0.1c 

8.72 

11770 

0.15 

12.74 

12070 

0.22 

18.23 

12370 

0.30 

24.25 

12520 

0.45 

35.94 

| Strut  No.  6—  j 

*10560 

0.01 

0.95 

*12070 

0.0k 

3.31 

12370 

0.06 

4.85 

12670 

0.10 

7.89 

12970 

0.15 

11.57 

13270 

0.25 

18.84 

13430 

0.34 

25.32 

13580 

0.?4 

54.49 

g 

*  These  points  were  not  considered  amenable  to  the  linear  plot  try  Southwell 
because  of  the  lack  of  accuracy  at  the  lower  values  of  load  and  deflection 
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TABLE  II.  COMPARISON  OF  THE  THEORETICAL  CRITICAL  LOADS  WITH  THE  SLOPE 
VALUE  FROM  TEE  SOUTHWELL  PLOTS  (Ref.  6) 


c  deduced  from 
best-fitting 
line  in 
Figure  9 
(mm) 

P^  estimated 
from  slope  of 
best-fitting 
line  in  Fig.  9 
(kg) 

P  as  given 
by  theoretical 
formula 
(kg) 

0.005 

3712 

3790 

0.005 

5^53 

5^75 

0.005 

8590 

8645 

0.005 

8758 

8610 

0.003 

11220 

10980 

0.030 

HO9O 

10920 

0.010 

12815 

12780 

0.010 

13750 

!  13980 

(5) 


Estimated  value. 


Theoretical  value 
(Col.  4) 


.980 
•995 
0.994 
1.017 
1.022 
1.015 
1.003 
0.984 


where 

P  and  y  =  the  load  and  the  corresponding  midpoint  deflection, 
respectively 

P^  and  y^  =  initial  values  of  P  and  y,  respectively 

P  . ,  =  the  first  Euler  load 
cnt 

TTX 

a  sin  — -  =  a  constant  for  a  given  reference  noint  x. 

1  L 

The  straight  line  obtained  by  plotting 


as  an  ordinate  against  y  -  y^  as  an  abscissa  cuts  the  horizontal  axis  at 
the  distance 

tv*  rj 

from  the  origin,  and  the  inverse  slope  of  the  line  is  Pcrit  -  P^. 

With  this  formula  he  was  able  to  analyze  some  of  the  von  Karmen  results 
which  were  not  amenable  by  the  original  approach. 

This  process  of  analyzing  experimental  data  was  extended  by  Fisher10' 11 
to  deal  with  raoie  complex  conditions.  The  case  which  he  considered  is 
depicted  in  Figure  11.  In  this  diagram  the  line  AA'  represents  the  center 
line  of  the  beam.  Thus,  the  conditions  assumed  for  the  purpose  of  analysis 
are: 

1.  Compressive  loads  P  are  applied  at  fixed  points  in  the  end  cross 
sections  of  the  beam. 

2.  The  end  eccentricity  of  P  is  small. 

3*  The  center  line  of  the  beam  in  the  unstrained  case  is  represented 
by  y  =  n(x) . 

4.  Forces  proportional  to  P  are  applied  to  the  beam  in  the  direction 
0^,  and  end  couples  proportional  to  P  are  applied  in  the  plane  0^. 

5.  The  forces  described  in  (4)  above  and  their  associated  reactions 
are  defined  by  saying  that  they  produce  in  the  absence  of  other 
loads  c.  bending  moment  P  •  N(x)  at  a  distance  x  from  the  origin. 

6.  Constant  forces  parallel  to  Oy  are  applied  at  various  points  along 

the  beam,  and  constant  end  couples  are  applied  in  the  plane  0  . 

7.  The  system  described  in  (6)  above  is  defined  by  saying  that  these 
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forces,  together  with  their  appropriate  reactions,  produce 
a  bending  moment  M(x)  at  position  x  and  end  rotation  m(x). 

With  these  assumptions  the  equilibrium  equation  may  be  written  as: 

d2 

^  ^2x  {o(x)  +  m(x)|  -  P  {n(x)  +  m(x)  +  6(x)j  +  M(x)  +  P*N(x)  (l6) 

or 

|d2  +  p2J  6(x)  =  -u2  {n(x)  +  N(x)  +  m(x)|  =  -^2H(x)say  (17) 

where 

6(x)  is  the  added  deflection  due  to  P  and  its  associated  loads  and 

2  P 
^  ~  El* 

The  deflection  6(x)  is  determined  by  equation  (17)  and  by  the  condition 
6  =  0  at  x  ±  L.  It  is  to  be  noted  that  there  is  an  exact  correspondence 
among  y,  N,  and  m,  in  that  it  is  not  these  individual  values  which  affect 
6  but  only  the  sum  H.  Thus,  we  have  the  most  important  results  -  namely, 
initial  bowing,  end-eccentricity,  transverse  loads,  and  end  couples 
whether  proportional  to  P  or  constant,  all  combine  to  form  one  total  equiv¬ 
alent  eccentricity. 

Fisher  then  proceeded  with  a  rigorous  analysis  and  demonstrated  that  South¬ 
well's  expression  should  be  written  as 

P1  v  -  6  =  H(0).  Fh(0,z)  (18) 

where 

H(6>)  =  the  sum  of  the  normal  forces  on  the  beam 

F..(0, z)  =  a  function  dependent  upon  the  normal  force  distribution 
n 

V  =6 

p* 

Wow  the  variation  of  F„  as  z  passes  from  0  to  1  is  a  measure  of  the 
deviation  from  linearity.  This  he  demonstrated  to  be  small  in  general. 

A  second  approximation  is  also  given  to  the  v/6  relationship.  This  is  in 
the  form  of  a  hyperbola. 

(v  -  5  -  Cl)(v  -  K  +  °l  -  0=  -°l(cl  -  Co)  <19> 


where 


S 

§ 

I 
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Figure  13 .  Plot  of  Lateral  Load  R  vs  Intercept  Value 
C  From  Figure  12. 
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P,  v,  and  6  have  the  significance  previously  ascribed,  and  £.  and  C 
are  constants.  K  has  the  value  0.123234.  ° 

To  test  his  theory,  Fisher  performed  a  series  of  experiments.  The  results 
are  portrayed  in  Figures  12,  13,  and  14.  It  is  seen  from  these  figures 
that  even  large  values  of  side  force  do  not  destroy  the  linearity  of  the 
6/p  versus  &  curves.  Thus,  the  tests  provide  substantial  verification  of 
Fisher's  analysis.  It  is  interesting  to  note,  too,  that  when  the  South- 
well  eccentricities  are  expressed  in  terms  of  the  side  force  values,  the 
relationship  is  linear  (Figure  13) .  This  is  clear  evidence  that  the 
intercept  with  the  axis  has  the  meaning  normally  ascribed  to  it. 

7 

In  his  paper  on  the  applicability  of  the  Southwell  plot,  Donnell  also 
considered  analytically  the  case  of  a  hinged  strut  with  continuous  elastic 
support.  This  strut  system  is  shown  in  Figure  15.  For  the  purpose  of 
his  analysis  he  assumed  an  initial  deformation  given  by 


CO 


n=l 


(20) 


and  that  the  movement  during  buckling  could  be  represented  by 


I  Wn  *  Sin  IT  (2l) 

n=l 

Now  the  total  energy  change  due  to  a  virtual  displacement  dW  must  vanish, 
or  n 


dW 

n 


j~EI/d2w\2 

L2\  2/ 
dx 


dx 


(22) 


=  P 


+w+dW 

n 


where 

El  =  the  bending  rigidity 
S  =  the  modulus  of  the  elastic  support. 
Thus,  it  can  be  shown  that 


Wn  n^EI  L2g  \ 

W'n  +  Wn  1  L2  n2*2^ 


W 


W 


n _ 

+  W 


n  n 


(23) 
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•fS**-  V*  '  v'  -r*x-  »-^r  *»r  w>v  • 


-*-  y  ^  *Jr-  ***-10  r ,  *  V^V*^  » 


This  equation  can  he  written  in  the  form 

V 


W  = 


n  P 


n 


n  -1 


(24) 


where 

Pn  =  the  classical  critical  loads  for  a  perfect  strut  having  a 
continuous  elastic  support. 

Once  again  the  relationship  between  the  elastic  deformation  and  the  applied 
load  reveals  a  form  similar  to  the  Southwell  representation  for  the  strut. 
Of  course,  in  a  problem  of  this  type,  which  may  easily  exhibit  modal  shapes 
having  harmonics  higher  than  the  fundamental  one,  it  is  important  to 
observe  and  record  adequately  the  buckling  shape  as  well  as  the  load  and 
deflection  parameters. 
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Hayashi  and  Kihira  in  a  report  recently  presented  to  the  Japan  Congress 
on  Testing  Materials  provide  some  experimental  results  for  this  problem. 

The  test  setup  which  they  used  is  depicted  in  Figure  16,  while  their 
results  are  portrayed  in  Figures  17,  18,  and  19 . 

Just  as  a  restraining  spring  located  along  the  length  of  a  column  causes 
an  increase  in  the  critical  load  reflected  by  a  change  in  slope  of  the 
Southwell  line  (as  the  above  referenced  tests  show),  so  does  a  point  of 
local  stiffness  reduction  caused  by  slight  yielding  of  the  material  in  a 
local  area  cause  a  change  in  the  apparent  critical  load.  This  situation 
is  clearly  seen  in  the  test  data  of  Fisher  given  in  Figure  20.  This  is  an 
important  facet  since  such  local  yielding  is  not  always  readily  apparent 
from  the  load  displacement  curve. 

The  discussion,  so  far,  has  centered  around  uniform  pin-ended  columns  and 
their  fundamental  critical  loads.  Such  restrictions,  however,  are  not 
necessary.  The  process  applies  for  all  end  conditions,  as  has  bepn  demon¬ 
strated  analytically  by  Ariaratnam^  and  experimentally  by  Hill.  4  The 
question  of  critical  loads  beyond  the  first  was  treated  by  Donnell 
This  analysis  was  perhaps  his  most  important  contribution  to  the  column 
problem. 

In  the  normal  Southwell  derivation,  as  in  Donnell's,  the  total  displacement 
at  any  point  is  the  algebraic  sum  of  the  various  harmonic  components  of  the 
deformation  at  that  point.  If  three  of  these  are  considered,  as  depicted 
in  Figure  21,  then  the  overall  deflection  6  at  the  1/3  and  2/3  positions 
is  related  to  the  harmonic  amplitudes  as  follows: 


6l/3  =  0.866(W1  +  w2) 
62/3  =  o.866(W1  -  W2) 
27 


(25) 

(26) 


Figure 


|P 


P 


16.  Experimental  Test  Setup  for  Column  Testing 
(from  Reference  12). 
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6  (mm) 


Figure  17.  Plots  of  P  vs  6  for  the  Buckling  of  Elastically- 
Supported  Columns. 
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Figure  l8.  Plots  of  6  vs  6/p  From  Data  Shown  in  Fig.  17 
(from  Reference  12). 
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in  which  expression  there  is  an  error  of  W_  and  higher  harmonics;  whereas  in 
the  second  approximation  formula  given  above,  equation  (28),  the  error  is 
and  higher  harmonics. 

It  is  clear  that  in  cases  in  which  a  large  number  of  harmonics  are  concerned, 
the  use  of  the  many  simultaneous  equations  would  add  greatly  to  the  data 
processing  as  well  as  to  the  acquisition.  But  by  using  the  principles  of 
Harmonic  Analysis,  it  can  be  shown  that  if 


^l/m’  62/W  ^(m  _  i/m) 


are  deflections  at  points  dividing  the  length  into  m  aliquot  parts. 


Wn  T  m  L6(l/m)Sin  IT  +  6(2/m) 


sin 


2nrr 

m 


+  ****•' 


%  -  l/m)siQ  <“  -  X>  f]  <2« 

Donnell  remarked  that  it  was  unfortunate  that  no  experimental  evidence 
existed  on  this  question.  However,  in  a  report  published  in  the  Journal  of 
Research  of  the  National  Bureau  of  Standards,  Tuckerman1^  examined  the 
Southwell  plot  from  the  point  of  view  of  the  Westergaard-^  general  theory, 
and  produced  data  to  verify  the  harmonic  theory.  This  experimental  work  was 
performed  by  McPherson  and  Levy.  From  this  laboratory  study,  he  was  able 
to  determine  first  and  second  critical  loads  for  one  column  tested  and 
first  and  third  for  another. 

The  experimental  results  are  presented  in  Table  III,  and  the  final  curves 
are  given  in  Figure  22  for  the  first  and  second  mode  problem.  The  theor¬ 
etical  critical  load  levels  are  compared  with  the  experimentally  determined 
values  in  Table  IV. 

The  corresponding  data  with  respect  to  the  determination  of  the  first  and 
third  critical  levels  are  given  in  Tables  V  and  VI  and  in  Figures  23  and  2k. 

It  should  be  noted  that  in  each  case,  the  agreement  between  the  first 
critical  load  and  the  Euler  load  is  good  (3  percent),  but  the  e.greement 
between  the  second  and  third  values  and  their  corresponding  theoretical 
counterparts  is  poorer  (ll  percent  and  12  percent). 
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TABLE  III.  CRITICAL  LOAD  OF  13-5-BY-1/4-EY-1/2-INCH  COLD-ROLLED  STEEL 
_ COLUMN  WITH  EQUAL  AND  OPPOSITE  ECCENTRICITIES  AT  THE  ENDS* 


c - 

. . . — - -  -  -  "  "  1 

Load 

AT** 

AM** 

AB** 

P1 

2a*** 

id  *** 
2°1 

P1 

d  *** 

pi 

lb 

io_lf 

xo'k 

10'^ 

lb 

10"^ 

10"1* 

lb 

10-1*  lo'* 

lb 

mm 

0 

0 

0 

0 

0 

- 

0 

-1.02 

-0.73 

0.01 

-40 

-0.37 

.00912 

-0.51  .01288 

■SB 

-1.93 

-1.20 

0.13 

-80 

-0.60 

.00750 

-1.03  .01288 

670— 

-2.73 

-1.63 

0.33 

-120 

-0.81 

.00679 

-1.53  .01275 

630— 

-3.40 

-1.93 

0.54 

-160 

-0.96 

.00602 

-1.97  .01231 

590- 

-4.10 

-2.19 

0.91 

-200 

-1.09 

.00548 

-2.50  .01252 

550- 

-4.67 

-2.34 

1.21 

-240 

-1.17 

.00488 

-2.94  .01225 

510— 

-5.23 

-2.53 

1.52 

-280 

-1.26 

.00452 

-3.38  .01205 

470— 

-5.73 

-2.66 

1.82 

-320 

-1.33 

.oo4i6 

-3.78  .01179 

430— 

-6.29 

-2.78 

2.19 

-360 

-1.39 

.00386 

-4.24  .01178 

390- 

-6.74 

-2.87 

2.52 

-400 

-1.44 

.00359 

-4.63  .01158 

350- 

-7.22 

-3.01 

2.83 

-440 

-1.50 

.00342 

-5 .02  .01142 

310- 

-7.69 

-3.05 

3-14 

-480 

-1.52 

.00318 

-5.42  .01130 

ro 

3 

1 

1 

-8.14 

-3.13 

3.47 

-520 

-1.56 

.00301 

-5.80  .01115 

230— 

-8.53 

-3.19 

3-79 

-560 

-1.60 

.00285 

-6.16  .01100 

£ 

O 

1 

1 

-9.01 

-3.26 

4.08 

-600 

-1.63 

.00272 

-6.55  .01090 

150— 

-9.40 

-3.35 

4.40 

-640 

-1.68 

.00262 

-6.9c  .01080 

110— 

-9.78 

-3.36 

4.71 

-680 

-1.68 

.00247 

-7.24  .01065 

*  This  table  reproduced  from  Reference  15. 

**  A  is  the  difference  in  strain  on  opposite  sides  of  the  specimen 

(bending  strain)  arranged  to  be  zero  at  the  arbitrary  "zero  load" 
for  computations  of  790  lb. 

***  d^  and  d?  are  the  bending  strains  corresponding  to  the  first  and 
second  critical  loads,  respectively. 

In 
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/AM\ 
L  “  \  2/ 

,'AT  -  AB^ 

Q2  "  V  2  ) 

. 

r——  | 
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Figure  22.  Lundquist-Type  Plots  for  First  and  Second  Modes 

of  <~n  Axially  Compressed  Column  (from  Reference  15 ). 
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FIRST  CRITICAL  LOAD 


From  strain  at  middle . 994  It 

From  strain  at  quarter  points . 1,015  lb 

Tl^EI 

rm,  _ _ "  A  A  i  AO  A  tv 


Theoretical 


1,020  lb 


£  P 

(Assumed  E  -  29  x  10  lb/in.  ) 


SECOND  CRITICAL  LOAD 

From  strain  at  quarter  points  . 

Theoretical  —■ . . 

.  .  .  3,700  lb 

.  .  .  4,100  lb 

L  /■  p 

(Assumed  E  -  29  x  10  lb/in.  ) 

*  This  table  reproduced  from  Reference  15 


TABLE  V.  CRITICAL  LOAD  OF  27.7-BY-3/o-BY~3/4-IN0H  COLD-ROLLED  STEEL  COLUMN  WITH  EQUAL 
ECCENTRICITIES  AT  THE  ENDS* _ _ 
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The  next  application  of  interest  appears  in  the  study  reported  hy  Bridget, 
Jerome,  and  Vosseller.1?  This  work  was  concerned  with  the  stability  of 
angle  struts.  The  instability  mode  for  such  structures  is  markedly 
dependent  on  the  width  of  the  flanges.  When  the  flanges  are  narrow,  the 
mode  is  of  the  Euler  Column  type,  but  when  they  are  wide,  local  or  torsional 
instability  can  occur.  It  was  found  in  the  investigation  that  the  South- 
well  plot  was  applicable  irrespective  of  mode.  The  overall  results  of  the 
study  are  given  in  Figure  25 .  A  typical  Southwell  plot  is  shown  in  Figure 
26. 

7 

These  experiments  were  conducted  under  the  guidance  of  Donnell,  who  pro¬ 
vided  theoretical  justification  for  the  procedure  in  a  report  published 
in  1938-  He  analyzed  the  stability  of  a  compressed  flat  panel  hinged  on 
three  sides  and  free  on  the  fourth,  as  shown  in  Figure  27.  Because  the 
deformed  shape  of  such  a  body  is,  to  a  fair  approximation,  a  developable 
surface,  extensional  stresses  can  be  neglected  in  the  internal  energy, 
and  the  initial  deformation  of  the  panel  can  be  represented  by 


w' 


s 

b 


l 


sin 


nrrx 

L 


and  the  additional  deflection  can  be  represented  by 


v  =  r-  )  W  sin  nnx 
b  L  n  -v— 


(30) 


Then  the  condition  where  the  energy  change  due  to  a  virtual  displacement 

3 


dW  vanishes  is 
n 


dW 


Et“ 


b  „  L 


n  aw-  ^  -r  J  i  (V  +  2(1  -  *>  (&>‘ 

n  12(1  -  p.  )  o  0  ox  os 


,2  2 


-.2  J2  n 
d  w  d  w] 

+  2* 


p  r*  ^  |»  L 

dx  ds  =  £  ! 

b  ^  J 
o  0 


w1  +  w  +  dW  ^  sin 

n  b  L  /J 


dx 


[(W  +  v)j  } 


dx  ds 


(31) 


where  E,  n,  and  t  are  the  elastic  modulus,  Poisson's  ratio,  and  thickness, 
respectively. 


Substituting  equation  (30)  into  (31),  we  obtain 
W-  r  Et3b  1  /n2!!2  6(1  - 


n  [  EtJb  1  /n_D_ 

W  '  +  W  LoM  2J  \  2 

n  n  12(1  -  [x  )  L 


w  W 
a _ 5L 

/  W  '  + 


W  Pn 
n  n 


hi 


P 


(32) 


"fch. 

where  P  is  the  n  critical  load.  This  expression  if  clearly  identical 
in  form  to  that  derived  by  Southwell®  and  suggested  by  Ayrton  and  Perry. 5 

If  the  linear  plot  method  (Southwell  Plot)  described  had  significance 
merely  with  regard  to  the  strut,  then,  without  question,  it  would  merit 
discussion.  But  the  value  of  the  technique  is  much  wider.  We  shall  deal 
with  various  successful  applications  of  the  procedure  in  chronological 
order. 
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In  1933  Gough  and  Cox'10  made  experiments  on  the  shear  buckling  of  plates 
with  a  view  to  checking  the  analytical  studies  of  Southwell  and  Skan.1^ 

In  the  first  series  of  tests  which  they  made,  the  buckling  loads  were 
determined  by  distortion  of  images  reflected  in  the  surface  of  the  strip. 
The  results  of  this  approach  were  very  positive.  In  all  cases,  there  was 
a  serious  discrepancy  between  theory  and  experiment.  The  "critical  loads" 
determined  as  described  were  found  to  be  representable  by  an  expression 
which  differed  from  the  theoretical  not  only  in  magnitude  but  also  in  form. 

As  a  consequence,  the  interpretation  of  the  tests  was  made  on  the  basis  of 
the  Southwell”  argument,  which  proceeds  as  follows. 

If  initial  irregularities  are  present  in  the  plane  of  the  sheet,  there 
will  be  a  component  of  the  wave  form  corresponding  to  the  first  critical 
load,  as  defined  by  theory.  If  the  amplitude  of  this  component  is  W  ,  the 
effect  of  a  shear  load  S  is  to  multiply  this  amplitude  in  the  ratio 
S/(S  -  S),  S  being  the  critical  load.  Thus,  the  corresponding  values  of 
S  anS  W  (elastic  deflection)  should  satisfy  the  relationship. 

W  =  SWq  /  (SQ  -  S)  (33) 

,-S 

i.e.,  ’  W  Jj-jjr  -  1 j  =  constant  (3^) 

an  equation  of  precisely  the  same  form  as  equation  (14), 

Thus,  the  plot  of  W  versus  W/S  should  be  linear,  and  the  slope  should 
correspond  to  the  first  critical  load  for  a  perfect  plate,  since  the  inter¬ 
cept  is  clearly  the  value  of  the  effective  imperfection. 

The  results  of  this  approach  are  portrayed  in  Figure  28.  Good  linear 
relationships  are  obtained.  The  conparison  of  the  theoretical  loads  and 
the  values  estimated  from  experiment  is  given  in  Table  VII. 

The  application  of  the  same  techniques  for  the  interpretation  of  the  2q 
instability  results  of  flat  plates  is  also  feasible.  Timoshenko  and  Gere 
determined  a  linear  relationship  similar  to  equation  (14)  for  plates  in 
which  there  is  no  midplane  stretching.  The  initial  deflection  must  be 
small  in  conparison  to  the  plate  thickness.  The  fourth-order  differential 
equation  of  equilibrium  is 


Ratio  of  wave  depth  to  shear  load 

Strip  DGLI,  150  -  l-l/2-in.  width  (w/s  plotted  from  0  as  zero). 
Strip  DGLI,  157  -  l-l/4-in.  width  (v/s  plotted  from  10  as  zero). 
Strip  DGLI,  156  -  1-in.  width  (w/s  plotted  from  10  as  zero). 
Strip  DGLI,  155  -  3A“i&*  width  (w/s  plotted  from  20  as  zero). 
Strip  DGLI,  15 4  -  l/2-in.  width  (w/s  plotted  from  30  as  zero). 

Figure  28,  Determination  of  Critical  Load  From  Gough  and  Cox 

Data  on  Thin  Strips  Under  Shear  (from  Reference  18). 
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TABLE  VII .  DETAILS  OF  STRIPS  TESTED  AND  MEASURED  VALUE!  OF  WAVE  DEPTH* 

Reference  mark 
of  strip 

DGLI  154 

DGLI  155 

DGLI  156 

DGLI  157 

DGLI  150 

Width  (in.) 

1/2 

3A 

1 

1  1/4 

1  1/2 

Length  (in.) 

10 

15 

20 

25 

30 

Thickness  (in.) 

0.0126 

0.0126 

0.0126 

0.0126 

0.0126 

Theoretical  buck¬ 
ling  load  (tons) 

16.3 

10.9 

8.2 

6.5 

5-4 

Buckling  load 
estimated  from 
wave -depth  read¬ 
ings  (tons) 

- 

11.4 

6.65 

6.86 

5.29 

s  estimated 

e  buckling  load 

1.045 

.81 

1.055 

.98 

S  ”  theoretical 
°  buckling  load 

Shear  load  (tons) 

Wave-Depth  (deflection  meter  reading)* 

0.2 

0 

0 

0 

0 

- 

0.25 

- 

- 

- 

- 

0 

0.4 

0 

6 

2 

1 

- 

0.5 

- 

- 

- 

- 

1 

0.6 

1 

7 

5 

2 

- 

0.75 

- 

“ 

- 

6 

0.8 

2 

o 

8 

3 

- 

1.0 

3 

12 

10 

4 

10 

1.2 

6 

14 

11 

6 

- 

1.25 

- 

“ 

- 

- 

12 

1.4 

8 

15 

14 

9 

- 

1.5 

- 

- 

- 

- 

18 

1.6 

10 

18 

18 

11 

- 

1.75 

- 

- 

- 

- 

20 

1.8 


2.0 


12 

16 


20 

21 


21 

24 


12 

14 


25 


TABLE  VII  -  Continued 


2.2 

20 

23 

29 

17 

•m 

2.25 

- 

- 

- 

- 

31 

2.4 

23 

26 

33 

19 

- 

2.5 

- 

- 

- 

- 

39 

2.6 

27 

29 

40 

22 

- 

2.75 

- 

' 

- 

- 

48 

2.8 

30 

31 

47 

25 

- 

3.0 

35 

35 

53 

27 

58 

3.2 

39 

39 

62 

32 

- 

3.25 

- 

- 

- 

- 

70 

3.4 

45 

42 

74 

35 

- 

3-5 

- 

- 

- 

- 

90 

3.6 

50 

48 

89 

40 

- 

3-75 

- 

- 

- 

115 

3-8 

57 

52 

104 

45 

- 

4.0 

65 

60 

127 

52 

150 

4.25 

- 

- 

- 

- 

195 

4.5 

- 

- 

- 

- 

252 

4.75 

- 

- 

- 

- 

323 

5.0 

- 

- 

- 

- 

401 

*  Data  repooduced  from  Gough  and  Cox  Reference  18. 

**  The  values  of  deflection  readings  used  for  determining  the  critical 
loads  are  those  included  between  the  horizontal  lines  marked  in  the 
table . 
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dx 
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ay 


(35) 


where 


w  =  elastic  lateral  displacement  of  the  median  surface 

v  =  initial  lateral  deformation  of  the  median  surface 
o 

D  =  flexural  stiffness 

q  =  transverse  loading 

and  the  in-plane  stresses  are 

a  -  N 
x  x 


CT  =  N 

y  jl 

t 

CT  =  N 
xy  xy 

t 

being  positive  in  compression.  From  the  equation  above,  it  was  shown  by 
Timoshenko  and  Gere  that  an  initial  curvature  produces  a  deflection 
identical  to  that  given  by  a  lateral  load  of  intensity. 


(36) 


~2 
9  w 


2 

o  w 


q  =  N 
o  x 


9x 


+  2N 
2  xy 


+  N_ 


dxdy 


.2 

o  w 


Sy 


(37) 


This  mathematical  identity  is  true  only  as  long  as  w  is  so  small  that  the 
principle  of  superposition  is  not  invalidated.  This°result  is  a  particular 
example  of  the  general  case  derived  by  Fisher  in  Reference  10.  (See  also 
Figure  36.) 

For  this  linear  plate  theory  under  the  restriction  of  a  uniaxial  load  with 
simply  supported  edges,  Figure  29,  the  elastic  lateral  deflection  may  be 
taken  as 


00  CO 


-I  1 


-  .  mTfx  .  nTTy 

B  sin -  sm  - - -*■ 

mn  a  b 


(38) 


m=l  n=l 

and  then  the  initial  deflection  is  defined  as 

k9 


1 


- 


ZXrf*  ■  '•*  ‘*r?^^'t*>*C*zv**  /»■'**$>■»»  ^*5*  « 


\-l  IV^T^? 


m^l  i=l 


Substitution  of  these  Fourier  relationships  into  equation  (37)  results  in 


B  \(z)\  (£)2]  **  .  ±j »  B  »)*♦  n  a  (a 

irn  L\a/  \b/  J  DL  x  mn\  a  /  x  mn\  a 


Thus, 


A  N 
mn  x 


4  [»2  *  &ff-  \ 


m  =  1,  2, 

n  =  1,  2, 


For  a  perfect  flat  plate  with  ideal  boundaries,  the  critical  uniaxial  loads 
per  unit  length  are  veil  established.  They  are 


W  =  N 
x  x 


cr  1  m,n  m  a 


=  4  [ma  +  (-f| 


in  =  1,  2,  . . . . 
n  =  1,  2,  . . . , 


Therefore,  equation  ( 4l)  becomes 


mn  "  N’ 


This  can  also  be  revrit.ten  as 


Bmn  L  N 


1  1  =  A 
j  mn 


m  =  1,  2,  3, 
n  =  1,  2,  3; 


JzlStliZU  1 


Clearly  this  result  is  in  a  one-to-one  correspondence  vith  that  established 
hy  Southwell  for  the  case  of  the  column. 

If  we  consider  a  realistic  square  plate  under  normal  test  conditions,  all 
experimental  evidence  indicates  that  the  deflection  is  of  the  first  modal 
shape  at  the  commencement  of  buckling,  and  equations  (42)  and  (44)  become 


(45) 


Equation  (46)  is  identical  to  equation  (l4).  Therefore,  the  Southwell 
Plot  is  also  applicable  in  this  case. 

As  buckling  develops,  however,  modal  changes  may  occur,  and  then 
consideration  must  be  given  to  the  effects  of  further  terms  which  result 
from  the  higher  harmonics.  Thus,  for  interpretation  of  test  data  appropr¬ 
iate  to  these  cases,  it  is  clearly  necessary  for  the  observer  to  describe 
and  record  accurately  the  modal  behavior  of  the  test  vehicle  and  to 
realize  the  need  for  careful  positioning  of  the  displacement  sensor.  In 
general,  for  a  square  plate  the  deformations  are  sinusoidal,  and  the  modes 
higher  than  the  first  do  not  normally  occur.  Hence,  in  practice,  the 
elastic  deflection  6  may  be  measured  at  almost  any  position  on  the  plate . 
Of  course,  for  maximum  sensitivity  the  center  point  should  be  chosen. 

7 

In  his  1938  paper,  Donnell  applied  the  finite  displacement  theory,  which 
he  had  previously  derived  in  1934^  for  thin-walled  cylindrical  shells,  to 
a  panel,  and  so  permitted  edge  restraint  and  allowed  for  in-plane  stretch¬ 
ing  of  the  plate.  On  this  basis,  if  $  is  the  usually  defined  stress 
function,  he  found  that 


1  A.  1  92»-  ^  r 
S,,  =  V*K 

9x 


r/  a2  \2 

Md)  - 


2  2  . 

S  w  5  w  1 


where  r  is  the  radius  of  curvature  and  x  and  s  are  longitudinal  and 
circumferential  coordinates.  If  $  and  w  are  harmonic  functions  of  x 
and  s  ,  the  internal  strain  energy  is 


rr  1 

v  =  i  r 


JJ  2LL2(1  - 


^2w) +  fe(v2$)  ]dxds 


In  the  case  of  a  flat  panel,  the  added  raidplane  extensional  terras  are 
those  of  the  first  order  squared.  These  are  due  to  the  elemental  rotations 


as  shown  in  Reference  21.  Let  the  initial  curvature  and  the  elastic 
movement  of  the  square  "flat"  panel,  hinged  on  four  sides  and  shown 
in  Figure  30,  be  respectively 


,  ...  .  TTx  .  TTx 

w'  =  W'sin  -=—  sin  — 
L  L 


...  TTx  .  TTx 
w  =  Tv  sin  —  sin  — 
L 


(49) 


Setting  l/r  =  0  shows  that  equation  (  47)  is  similar  to  the  compatibility 
equation  for  the  flat  plate  "large  displacement"  theory  derived  in  1910 
by  von  Karman.22  Now  further  substituting 


K  = 


1  +  2 


and  equation  (49) 


into  equation  (46)  results  in 


= 


Et 

L 


TT^ 


/,  _W'\  .2  (  2  TTx 

^1  +  2— J  w  (^cos  —  cos 


2  Its 


-  sin 


TTx  .2  TTsN 
L  S1  L/ 


EtTT^ 


+  2W 


)W  ( 


2TTx 

COS  -7—  +  cos 
Jj 


^6) 


(50) 


When  possible  conditions  for  displacements  in  the  plane  of  the  plate  at  the 
edges  are  neglected,  equilibrium  in  the  plane  of  the  plate  is  evidently 
satisfied  if 


’2$  =  (W  +  2W')w(cos  ~  +  cos  ~) 

Sir  \  L  L  / 


(5l) 


The  condition  that  the  total  work  due  to  a  virtual  displacement  dW  vanishes 

o  o 


“  Lai  (w'  +  v)]  }  dxds 


(52) 


Substituting  equations  (48),  (49)  ,  and  (51)  into  equation  (52)  results 
in 


P  = 


W 

W'  +  w 


[-  ■  — ~ - ][l  +  (W  +  2W')(W  +  W')J 

3(1  -  n  )  L  8t 


53 


=  w  +  w  pc  L1  + 


-H-i  (W  +  2W*)(W  +V)]  (53) 


where  P  is  the  critical  load  without  initial  curvature, 
c 

It  is  clear  from  equation  (53)  that  if  W*  and  W  are  very  small  in  relation 
to  the  plate  thickness,  then  the  expression 

2 

f  =  yj.  +  2W»)(W  +  W;)]  (5^) 

8t 

is  negligible  and  for  all  practical  purposes  the  Southwell  approximation  is 
valid.  However,  as  the  ratios  of  W'/t  and/or  W/t  increase,  the  Southwell 
condition  is  violated.  Nevertheless,  if  the  ratio  of  W'/W  is  small, 
equation  (53)  can  be  written  as 


r  =  Pc[i-GB  +  ^(1  +  f.)] 


where 


which  may  be  further  approximated,  with  little  error,  to  the  expression 


P  =  P 


[l  +  yW2] 


Thus,  if  a  new  variable  is  chosen  such  as  Z  =  W  ,  equation  (56)  may  be 
written  as 


P  =  Pc(l  +  yz) 


which,  of  course,  is  a  linear  equation  in  the  variables  P  and  z.  This 
result  is  similar  in  form  to  that  published  by  Yoshiki  et  al.23>24,25 
and  also  bears  a  remarkable  similarity  to  Qusinec's^^  results  which  are 
discussed  later  in  this  report.  It  can  be  visualized  easily  that  i-r  a 
graphical  representation  of  equation  (57)  is  made  for  the  variables  P  and  z, 
then  the  straight  line  will  intercept  the  load  axis  at  a  value  of  P  which 
corresponds  to  the  theoretical  critical  load  value  Pc . 

At  this  stage  it  is  interesting  to  determine  the  bounds  of  applicability  of 
the  Southwell  approximation  in  terras  of  the  amplitude  of  permitted  motion. 
For  this,  the  important  term  in  Donnell's  equation,  equation  (53),  is  of 
course  the  expression  for  7,  equation  (54).  Thus,  if  a  valid  Southwell 
condition  is  to  exist,  this  coefficient  Y  must  be  very  small  in  comparison 


55 


to  unity.  In  general,  the  upper  hound  of  applicability  for  equation  (53) 
would  appear  to  occur  when  the  square  of  the  measured  displacement  is 


8t2 

3(1  -  n2) 


(58) 


where  Y  «<  1. 

Equation  (53)  can  be  rewritten  similar  to  the  Southwell  form  as 


p_  rw  +  wi 
p  L  w  J 

c 


1  +  Y 


or 

r(l  +  Y)PC  _  w  _  w,  (59) 

P  J 

Hence,  it  is  obvious  that  the  numerical  value  of  Y  is  the  measure  of  the 
precision  with  which  the  Southwell  Plot  will  predict  the  classical  critical 
load.  In  other  words,  the  coefficient  Y  is  the  error  introduced  into  the 
evaluation  if  no  bounds,  other  than  the  limitations  of  the  original  "large 
displacement"  derivation,  are  placed  on  W  (or  W/t) .  For  example,  if  the 
error  is  to  be  limited  to  1  percent,  then  Y  =  .01;  therefore,  the  upper 
bound  placed  on  W  is  given  from  equation  (58)  as 

(w)2  *  - ^-§-(.01)  (60) 

3(1  -  n2) 

Experiments  conducted  in  the  Aeronautical  Structures  Laboratory  at  Stanford 
University  have  demonstrated  the  validity  of  the  Southwell  approximation 
for  plates  with  restrained  loaded  edges  and  simply  supported  unloaded  edges. 
These  results  are  shown  in  Table  VIII  and  Table  IX  and  in  Figures  31,  32, 

33,  3^\>  and  35*  They  have  shown  that  imperfection  magnitude  and  side  rail 
clearance  are  linearly  related,  as  are  imperfection  values  and  the  intensity 
of  a  central  normal  force.  This  result  matches  that  of  Fisher^  for  a 
column. 

If  further  justification  of  the  Southwell  linear  approach  for  flat  plates 
is  needed,  it  can  be  obtained  by  reevaluation  of  the  analytical  work  of 
Coan^°  and  the  experimental  results  of  Hoff,  Boley,  and  Coan.^7  Coan 
analytically  investigated  the  buckling  of  rectangular  isotropic  plates 
having  small  initial  curvature,  simply  supported  loaded  edges,  and  simply 
supported  unloaded  edges  which  were  free  to  warp.  The  initial  curvature 
was  restricted  to  be  small.  The  theoretical  curve  relating 


P  /P 
x' 


cr 
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JTABLS  VIII.  BEHAVIOR  OF  A  COMPRESSED  4-PLY  FIBERGLAS  PANEL  OF  181  CLOTH  ■ 
6-IN.-BY-6-IN.  LOADED  EDGES  CLAMPED  AND  UNLOADED  EDGES  SIMPLY 
SUPPORTED  WITH  VARIOUS  SIDE  RAIL  CLEARANCES.  _ 


Load  (lb  ) 

Strain  p  in. /in. 

20 

20 

20 

20 

30 

23 

23 

55 

40 

30 

30 

67 

50 

38 

38 

85 

60 

48 

48 

108 

70 

62 

70 

137 

80 

83 

102 

185 

90 

115 

140 

250 

100 

158 

192 

339 

110 

219 

278 

445 

120 

300 

393 

570 

130 

400 

528 

67  6 

140 

508 

648 

775 

150 

615 

748 

892 

160 

720 

853 

960 

170 

818 

958 

1092 

180 

900 

1054 

1175 

Side  Edge  Clearance 

(in.) 

_ 

0.0015 

0.002 

0.004 

Tests  conducted  at  Aeronautical  Structures  Laboratory,,  Stanford 
University. 
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TABLE  IX.  BEHAVIOR  OF  A  COMPRESSED  4-PLY  FIBERGLAS  PANEL  OF  181  CLOTH  - 
6-IN.-BY-6-IN.  LOADED  EDGES  CLAMPED  AND  UNLOADED  EDGES  SIMPLY 
SUPPORTED  WITH  VARIOUS  NORMAL  FORCES  APPLIED  AT  CENTER  OF 

PANEL 

Load  (lb) 

Strain  Units 

0 

0 

0 

0 

0 

0 

20 

o.to 

0.80 

0.80 

0.80 

0.80 

to 

0.95 

i.to 

1.60 

1.80 

2.05 

60 

1.50 

2.15 

2.50 

3.10 

3.55 

80 

2.25 

3.25 

3.80 

4.85 

5.35 

100 

3.30 

4.60 

5.45 

7.20 

7-75 

120 

4.80 

6.55 

7.65 

10.15 

10.90 

ito 

6.95 

9.25 

10.80 

14.20 

15.00 

160 

10. to 

13.20 

15.30 

19.60 

20.30 

180 

16.20 

19.30 

21.40 

26. to 

26.90 

200 

25.  to 

27.40 

29.20 

3^30 

34.30 

Normal  Force 
(gm) 

Test  1 

0 

Test  2 

75 

Test  3 

130 

Test  4 

190 

Test  5 

250 

Tests  conducted  at  Aeronautical  Structures  Laboratory,  Stanford 
University. 


Buckling  Tests  of  Fiber  Glass  Panels  Conducted  at  Stanford  University  With 
Various  Side  Edge  Clearances. 
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From  Figures  32  and  3^ 


sod  Kq,  Figure  4  or  Reference  26,  where 
Px  =  average  compressive  load 
P  =  critical  load  for  a  perfect  plate 

K  =  elastic  deformation  at  the  plate  center  divided  by  the 
plate  thickness, 

was  used  to  obtain  the  data  listed  in  Table  X.  This  information  has  been 
plotted  in  the  Southwell  manner  and  is  shown  in  this  form  in  Figure  36. 
The  resulting  straight  line  has  an  inverse  slope  of 


P  =  jr-  =  o-96 


which  is  within  4  percent  of  the  classical  relationship. 

27 

In  the  experimental  work  of  Hoff,  Boley,  and  Coan,  "flat"  fiberglasr 
plates  were  buckled.  Lateral  deflections  and  strains  of  the  plates  were 
measured  and  load-deflection  curves  as  well  as  load-strain  curves  were 
shown.  Some  of  the  data  from  Figures  6,  9,  and  10  of  the  report  are 

shown  in  TabLes  XI  and  XII.  These  data  were  recast  in  the  Southwell  manner 
and  plotted  in  Figures  37,  38,  and  39- 

This  technique  can  be  applied  also  to  the  case  of  a  biaxially  loaded  plate. 
The  dimensions  and  coordinates  are  shown  in  Figure  40.  The  equilibrium 
equation  which  holds  true  under  the  assumption  of  a  Hookean  homogeneous, 
isotropic,  and  perfectly  flat  plate  is 

DvV  =  -  N  ii-S  £$  (62) 

X  C.  y  2 
°x  oy 

when  the  effect  of  an  initial  deformation,  w  ,  is  included,  this  equation 
becomes 

T>r>4  ,,  I  d^w  S^w  *j  .  fdjw  S^w  ”i  ,,  . 

DV  w  =  -  N  '  — *  +  _ o  -  N  — p  +  _ 0 J  (63) 

For  a  simply  supported  plate  where  the  biaxial  forces  N  and  N  are  uni¬ 
formly  applied  in  the  plane  of  the  middle  surface,  the  appropriate  elastic 
deflection  is 


B  sin 
mn 
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TABLE  X.  DEFLECTION  BEHAVIOR  OF  A  SIMPLY  SUPPORTED  SQUARE  PLATE  WITH 

SMALL  INITIAL  CURVATURE  LOADED  IN  EDGE  COMPRESSION,  UNIFORMLY 
DISPLACED  LOADING  EDGES  AND  STRESS-FREE  SUPPORTED  EDGES 

P 

K 

K 

_ c 

X 

— 

P 

cr 

c 

P 

X 

p 

cr 

0 

0 

0 

0.347 

.048 

0.1360 

0.533 

.072 

0.1680 

0.459 

.0848 

0.1850 

0.500 

.0983 

O.I966 

0.553 

.1200 

0.2170 

0.569 

.1280 

0.2250 

0.597 

.1380 

0.2310 

0.621 

.1600 

0.2580 

0.710 

.2240 

0.3150 

O.766 

.2720 

0.3550 

0.826 

•  3520 

0.4260 

O.87I 

.4320 

0.4960 

0.903 

.5040 

0.5580 

G.992 

.606 

0.6110 

Data  obtained  from  the ore 

tical  curve,  case  1,  Figur 

26 

e  4  of  Coan. 
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TABLE  XI.  LOAD  DISPLACEMENT  DATA  FOR  BUCKLING  TESTS  ON  SQUARE 

FIBERGLAS  PANELS 

^  3 

lb  x  10-5 

6 

in. 

6  -  6  /P  -  P 
o'  0 

in.  ---6 

IT*10 

0 

P4 

1 

P4 

6-6 

0 

Panel  (28a  +  he),  (Data  obtained  from  Fig.  S 

>,  Ref.  27) 

■SB 

.0167 

0 

0 

- 

.0222 

0.43 

.0055 

.0128 

-  :l 

.0278 

0.70 

.0111 

.01585 

4.31 

.0333 

0.95 

.0166 

.0175 

4.55 

.0431 

1.19 

.0250 

.0210 

4.78 

.0514 

1.42 

.0347 

.0245 

4.96 

■  .0625 

1.60 

.0458 

.0286 

5.21 

.0833 

1.85 

.0666 

.O36O 

5. to 

.1028 

2.04 

.0816 

.0422 

5-57 

.1223 

2.21 

.1056 

.0477 

5.76 

.1528 

2.40 

.1361 

.0567 

Panel  (24d/4c) 

,  (Data  obtained  from  Fig.  7> 

Ref.  27) 

4.52 

.00190 

0 

0 

- 

4.83 

.00240 

0.31 

.0005 

.00160 

5.31 

.00344 

0.79 

.00154 

.00195 

5.86 

.00495 

1.34 

.00305 

.00177 

6.38 

.00715 

1.86 

.00525 

.00282 

6.83 

.00990 

2.31 

.00800 

.00346 

6.92 

.01070 

2. to 

.00880 

.00366 

7.07 

.01237 

2.55 

.01047 

.00412 

7.32 

.01430 

2.80 

.01240 

.00443 

7.65 

.01730 

3.13 

.01540 

.00492 

8.07 

.02090 

3.55 

.01900 

.00535 

8.52 

.02450 

4.00 

.02260 

.00565 

Panel  (28a  +  4c),  snug  side  edge  clearance, 

(Data 

obtained  from 

Fig.  10,  Ref. 

27)  Test  one. 

5.00 

0 

0 

• 

5.24 

0.24 

.0026 

.0108 

5.49 

0.49 

.0066 

.0135 

5.64 

.0347 

0.64 

.0120 

.0187 

5.78 

.0373 

O.78 

.0146 

.0187 

5*90 

.0427 

0.90 

.0200 

.0222 

6.07 

.0520 

1.07 

.0293 

.0274 

6.20 

.0627 

1.20 

.otoo 

.0333 

6.35 

.0800 

1.35 

.0573 

.0424 

6.51 

.0973 

1.51 

.0746 

.0494 

6.68 

.1213 

1.68 

.0986 

.0587 

6.82 

.1427 

1.82 

.1200 

.0660 
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TABLE  XI  -  Continued 


68 


Thus,  the  initial  deformation  is  assumed  to  "be 


where  m  and  n  represent  the  number  of  buckles  or  the  number  of  half 
waves  in  the  x  and  y  coordinate  directions,  respectively.  Substituting 
these  Fourier  relationships  into  equation  ( 63)  leads  to 


for 


m  =  1,  2,  . 


n  =  1,  2,  . 

This  equation  ( 67),  relating  the  amplitude  of  the  elastic  deformation, 
the  amplitude  of  the  assumed  inperfection  (a  constant),  the  classical 
critical  loading,  and  the  biaxial  loads,  is  similar  in  type  to  equation 
( 4l)  for  the  uniaxial  loading  case. 


For  the  perfect  flat  plate  with  ideal  boundaries,  the  critical  biaxial 
loads  are  directly  dependent  upon  the  aspect  ratio  of  the  plate,  p  and  are 
determined  from  the  relation 
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(68) 


The  critical  combinations  for  N  and  N  for  any  particular  value  of  a/b  are 
determined  by  minimization  of  e^uation'y(68)  with  respect  to  m  and  n.  The 
usual  interaction  curve  clearly  outlining  the  stability  bounds  for  an  a/b  = 
1.0  is  shown  in  Figure  4l. 


If  a  square  plate  is  considered,  the  aspect  ratio  will  have  a  value  of 
unity,  a/b  =1.0.  Therefore,  equation  (68)  becomes 
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The  interaction  for  this  particular  case  having  the  nondimens ional  form 
of  coordinate  axes  of 
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is  shewn  in  Figure  4l.  From  this  graph,  the  classical  critical  loads  can 
be  determined  for  particular  biaxial  load  ratios.  Thus,  if  the  loading 
ratio 

N 

K  "  Ny 

is  substituted  into  equation  (69),  the  critical  loading  for  a  perfect  plate 
would  become 
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After  substitution  of  this  expression  into  equation  (67),  the  elastic 
deflection  is 
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(71) 


-  1 


where  N  =  N  KN  .  This  implies  that  the  loading  rate  in  the  two 
directions  will  maintain  a  constant  proportionality.  Transposing  equation 
(71)  results  in 
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(72) 


for  the  particular  m  and  n  which  correspond  to  Ncr  obtained  from  Figure  4l. 
If  the  values  of  the  total  in-plane  loads 

P  =  N  a 
x  x 

P  =  N  a 

y  y 

are  substituted  into  equation  (72),  we  obtain 
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Thus,  once  again  the  relationship  between  the  varieties  of  the  elastic 
deformation  and  the  load  shows  a  direct  correspondence  to  the  Southwell 
equation  {lb) .  The  experimental  work  to  check  this  derivation  should 
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be  completed  in  the  very  near  future. 


In  the  analysis  performed  above,  the  imperfection  of  the  plate  was  con¬ 
sidered  to  be  an  initial  deformation  whose  shape  was  similar  to  the  expect¬ 
ed  buckle  mode.  An  alternate  manner  in  which  to  consider  the  imperfection 
is  to  suppose  that  an  additional  transverse  load  is  applied  at  some  point 
on  the  plate.  This  technique  has  received  attention  by  Loo  and  Evan- 
Iwanowski^S  in  the  experimental  study  of  instability  of  spherical  caps. 

Their  work  is  discussed  in  our  paper  on  shells .29  it  is,  of  course, 
validated  by  the  tests  of  Fisher  and  by  the  tests  reported  here  on  panels. 

Turner,  Martin,  and  Weikel^,^‘  numerically  examined  a  biaxially  loaded 
plate  with  such  an  added  normal  force.  Their  computer  results  are  listed 
in  Table  XIII  and  depicted  graphically  in  Figure  42.  These  computed 
points  are  essentially  collinear.  This  result  tends  to  give  credence  to 
the  suggestion  that  a  knowledge  of  the  nature  of  the  basic  imperfection 
is  not  of  importance  if  it  is  desired  merely  to  confirm  the  theory  for 
perfect  bodies.  It  is  obviously  of  importance  in  the  compilation  of  data 
designed  to  help  in  the  formulation  of  practical  theories  of  real  struct¬ 
ures  . 

Real  structures,  of  course,  are  very  rarely  subjected  to  uniformly  dis¬ 
tributed  loadings  or  restrained  in  the  simple  manner  considered  so  far. 

Walker 32  has  made  a  substantial  effort  to  consider  flat  plates  under 
conditions  which  more  closely  approximate  those  encountered  in  actuality. 

His  study  has  been  directed  toward  the  establishment  of  a  stability  theory 
for  plates  axially  compressed  by  nonuniformly  distributed  load.  In  the 
analysis  performed,  a  wide  range  of  possible  boundary  conditions  wa3  taken 
into  account.  The  theoretical  treatment  was  evaluated  in  the  laboratory.  The 
tests  which  were  performed  and  reported  in  the  given  reference  appear  to 
be  very  thorough,  well  conducted  tests.  In  the  analysis  of  the  data  derived, 
Walker  demonstrated  that  a  linear  relationship  existed  between  6  and  6/p, 
but  he  offered  no  analysis  to  confirm  that  this  was  so.  Typical  results 
are  given  in  Figure  43.  In  presenting  the  results,  the  author  pointed  out 
that  the  correlation  of  the  theoretical  critical  load  value  and  the  insta¬ 
bility  level  predicted  from  the  analysis  of  the  test  is  very  good.  This 
can  be  seen  from  the  graphs  in  Figure  44. 

However,  it  must  be  admitted  that  while  flat  plates  do  play  an  important 
part  in  everyday  structures,  generally,  conditions  are  more  complex.  In 
aeronautics,  combinations  of  plate  and  column  (Usually  called  a  stringer 
in  this  application)  are  frequently  met.  The  testing  of  components  of  this 
type  is  fraught  with_all  the  difficulties  of  the  simpler  cases.  Ramberg, 
McPherson,  and  Levy*5"*  have  carried  out  an  experimental  study  of  this  problem. 
In  their  analysis  of  the  test  data,  they  applied  a  Southwell- type  plot.  They 
demonstrated,  as  is  clearly  seen  in  Figures  45  through  48,  that  it  gave  good 
results  for  stringers  attached  to  sheets,  irrespective  of  their  failing  mode- 
They  were  unable  to  apply  the  process  to  the  sheet  between  stringers  due  to 
the  lack  of  appropriate  data  prior  to  the  instability  value . 
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TABLE  nn.  COMPUTER  RESULTS  -  RECTANGULAR  PLATE  STABILITY  PROBLEM 


*  Reference  31 
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Southwell  Plot  for  Biaxially  Loaded  Plate.  Computed 
Data  Points  From  Reference  31 • 
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Figure  44.  Comparison  of  Experimental  Results  and  Theory.  Data 
Points  Depict  the  Critical  Load  Determined  From 
Southwell  Plots  From  Reference  32. 
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Figure  45.  Southwell  Plot  of  Twisting  Instability  of  Short 
Z-ltype  Stringer  From  Reference  33 • 
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Figure  46.  Southwell  Plot  of  Twisting  Instability  of  Stringer 
of  Specimen  6  From  Reference  33* 


Generally  speaking,  references  to  the  Southwell  Plot  define  the  displace¬ 
ment  parameter  as  a  motion  normal  to  the  direction  of  loading.  However,  it 
is  relatively  simple  to  demonstrate  that  there  is  a  one-to-one  correspondence 
between  this  deflection  and  some  strain  parameter  in  the  plane- of  the  body. 
The  first  report  to  utilize  this  point  was  that  of  Tuckerman.  ^  Indeed,  in 
the  harmonic  test  previously  discussed  with  respect  to  his  report,  the 
determination  of  the  Southwell  Plots  was  made  using  strains  on  the  surface 
of  the  column.  The  linear  relationship  between  normal  displacement  and 
elastic  bending  conponent  strains  for  flat  plates  is  demonstrated  here  in 
Figure  49  from  tests  made  in  conjunction  with  the  flat  plate  studies  con¬ 
ducted  at  Stanford  University  and  those  previously  mentioned.  The  use  of 
the  Southwell  Plots  based,  uppn-strain  measurements  also  has  been  outlined 
and  verified  by  Gregory  ^  3  5  ( for  application  to  struts. 

Thus,  concern  has  been  mainly  with  elastic  structural  instabilities  which 
have  been  produced  by  usual  force  systems .  However,  in  1961  Qneinec3 
analytically  studied  the  behavior  of  a  centrally  heated  "flat"  circular 
plate  with  free  edges  and  checked  these  results  experimentally.  In  the 
course  of  his  work,  he  treated  the  problem  using  both  the  linear  and 
nonlinear  theories.  A  most  important  step  was  taken  when  the  author 
correlated  the  test  results  with  theory.  He  revealed,  on  the  basis  of 
linear  plate  theory  which  included  the  possibility  of  a  constant  initial 
geometric  imperfection,  the  validity  of  the  Southwell  Plot  for  this  partic¬ 
ular  loading.  Also,  a  deflection-squared  law  (6^)  was  derived  using  the 
large  displacement  analysis  and  this  was  correlated  with  some  experimental 
data. 

The  basis  of  Queinec's  analyses  is  as  follows.  The  temperature  distributicn 
and  edge  conditions  were  assumed  to  be  perfectly  axisymmetric .  lne  plate 
was  considered  to  be  of  uniform  thickness,  at  uniform  temperature,  and 
free  from  stresses  in  its  initial  state.  The  temperature  variation  through 
the  thickness  was  neglected,  as  was  the  effect  of  gravity  forces.  The  usual 
bending  theory  assumptions  were  made. 

The  general  equations  of  plate  bending  (Karman  equations)  modified  for 
temperature  effects  were  used. 
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These  equations,  transformed  into  polar  coordinates  and  reduced  for  the 
axisymmetric  case,  are 
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(75) 


D  _2_2  1  d  faa  dv\ 

t  w  r  dr  \dr  dr) 


or  these  two  equations  can  he  written  as 
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the  usual  Airy  stress  function  for  circular  plates 
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Unfortunately,  these  coupled  nonlinear  differential  equations  have  no  closed 
form  solutions  except  for  very  special  cases;  i.e. 

(no  membrane  stresses)  or 

(plane  stress  problem) . 

So  Queinec  utilized  the  well-known  energy  method  for  determining  the 
approximate  critical  temperature  for  both  the  linear  and  nonlinear  theories. 

Because  his  experiments  revealed  that  lateral  deflections  occur  well  below 
the  critical  value  of  temperature,  he  realized  that  these  plates  were  not 
absolutely  flat.  Hence,  Queinec  considered  the  effect  of  initial  imper¬ 
fections.  The  initial  deflection  was  assumed  to  be  axlsymmetric  and  was 
defined  by  w^r).  When  the  initial  temperature  is  uniform  and  the  initial 
state  is  free  from  stresses,  the  additional  deflection  w(r)  is  a  solution 
of  the  differential  equation 
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where 

err  =  stress  due  to  temperature  rise  T(r)  =  T  (©) 

and  this  equation  is  valid  for  small  deflections.  Equation  (76)  is  no  longer 
homogeneous.  When  the  initial  imperfection  is  proportional  to  the  elastic 
deflection,  the  following  expression  is  obtained  from  equation  (77): 
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This  relationship  can  also  be  written  as 
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which  of  course  has  the  identical  graphical  form  of  the  Southwell  Plot. 

This  researcher  was  probably  the  first  to  investigate  the  case  when  the 
imperfection  and  the  additional  deflection  were  not  of  the  same  form.  Using 
the  energy  equation,  he  derived  the  following  expression: 


(79) 


where  p  =  r/a  (a  =  radius  of  plate) 


T  .  p  .dp 
T 


(80) 


The  subscripts  denote  particular  radial  positions;  thus,  T  is  the  temper¬ 
ature  at  r  =  0.  Queinec  then  showed  that  this  equation  reduces  to 
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where  the  coefficient  K,  depends  only  on  the  initial  shape.  Equation  (8l) 
reduces  once  again  to 

wo  Lf—  -  X]  -  h  •  “io  t82) 

O 

This  equation  then  is  similar  to  equations  (78)  and  (l^).  AH  three  are 
equal  to  constants. 

The  central  temperature  versus  the  central  deflection  or  specimens 
numbers  one  and  two  is  plotted  in  Figure  50.  The  corresponding  Southwell 
Plots  are  illustrated  in  Figure  51*  The  critical  temperature  for  the 
center  of  the  plate  determined  from  this  Southwell  Plot  is  T  =  59  F. 
This  shows  a  7.8-percent  difference  from  the  value  °crit 

determined  from  the  energy  approximation  which  is  normally  unconservative. 
Therefore,  this  percent  of  difference  should  actually  be  smaller  when 
compared  to  an  exact  solution. 

The  author  appears  to  have  been  the  first  to  realize  that  it  was  possible 
to  develop  a  relationship  between  the  postbuckling  deformations  and  the 
initial  buckling  load.  For  the  problem  which  he  considered,  he  derived  the 
following  equation: 
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where  T  and  w  are  the  midpoint  temperature  and  displacement,  respectively, 
Tcrit  ^  the  c^itical  temperature  and  Kg  is  a  constant. 

2jO 

In  a  report  published  in  1962  by  Fernandez-Sintes,  Horton,  and  Hoff, 
dealing  with  the  thermal  buckling  of  annular  plates,  the  Southwell  plot 
was  used  to  correlate  the  derived  theoretical  results  with  the  test  data. 

Moir£  fringe  techniques  were  used  to  measure  the  lateral  deflections  when 
the  plate  was  heated  in  an  axisymmetric  manner.  The  test  setup  is  shown 
in  Figures  52,  53;  and  5^- 

The  relationship  derived  by  Queinec  and  just  discussed,  equation  (82),  was 
utilized  by  these  researchers  to  compare  the  experimental  results.  Figures 
55,  56,  and  57  show  the  curves  of  central  edge  temperature  and  central  edge 
deflection  versus  time  for  the  three  test  specimens.  The  corresponding 
Southwell  plots  are  shown  in  Figures  58,  59;  and  60.  The  values  of  critical 
temperature  determined  from  these  graphs  were,  respectively, 
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The  comparative  theoretical  values  were,  respectively, 

<Wl  "  U9-5"P 

Merit’s  -  m’F 

<Tcrit>3  "  102‘P 

The  differences  of  these  results  were  all  within  the  order  of  accuracy 
predicted  for  the  particular  test  equipment. 

Once  again  the  value  of  the  Southwell  plot  for  correlation  of  a  theoxy 
has  been  shown.  The  accuracy  of  these  results  most  certainly  will  be 
improved  with  the  intended  implementation  of  noncontacting  measurement  dev¬ 
ices  of  high  sensitivity  and  accuracy. 

This  general  process  for  correlating  experiment  and  theory  can  also  be 
applied  to  the  case  of  the  lateral  instability  of  the  deep  beam  under  the 
action  of  a  concentrated  force  applied  in  the  vertical  plane.  Such  a  beam 
when  simply  supported  in  two  planes  is  shown  in  Figuies  6l  and  62.  In  this 
type  of  structure,  where  the  t/h  ratio  is  small,  beading  takes  place  in  the 
plane  of  the  greatest  flexural  rigidity,  and  lateral  buckling  may  occur  at 
a  certain  critical  value  of  the  applied  load.  Both  bending  and  twisting 
must  be  considered. 

The  Rayleigh-Ritz  method  is  applied  to  the  total  potential  energy  to  deter¬ 
mine  the  relationships  between  the  load  and  both  the  elastic  lateral  defor¬ 
mation  and  the  elastic  torsional  deformation  for  a  deep  beam  with  initial 
imperfections  simply  supported  in  two  planes. 

As  the  beam  buckles  laterally,  the  strain  energies  of  lateral  bending  and 
twist  about  the  longitudinal  axis  increase  while,  at  the  same  time,  a 
certain  amount  of  work  results  from  the  movement  of  the  point  of  application 
of  the  concentrated  load.  Assuming  that  the  small  change  in  the  strain 
energy  of  bending  in  the  plane  of  the  beam  during  buckling  is  negligible, 
we  find  that  the  change  in  the  strain  energy  due  to  the  lateral  buckling  is 
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Figure  60.  Deformed  Elemental  Cross  Section  of  Deep  Beam 

Simply  Supported  in  Two  Planes  With  a  Concentrated 
Load  at  the  Centroid. 
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polar  second  moment 

Poisson's  ratio 

the  lateral  deflection 

the  angle  of  rotation  of  the  cross  section. 


When  the  symmetry  of  the  buckled  shape  is  considered,  equation  (84)  may  be 
written  as 
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Now  the  work  done  by  the  applied  load  during  buckling  is 
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The  term 
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is  the  infinitely  small  arc  described  by  the  end  of  the  beam  with  respect 
to  the  elemental  cross  section.  It  Js  caused  by  the  bending  of  the  element 
in  the  T]§  plane.  The  vertical  component  of  this  arc  when  the  angle  of 
rotation,  tp,  is  small  is 
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This,  then,  is  the  small  distance  through  which  the  applied  force  moves. 

The  total  work  done  is  one-half  the  summation  across  the  beam  of  these  ele¬ 
mental  vertical  components  multiplied  by  the  applied  load.  When  symmetry 
of  the  deformed  shape  is  assumed,  this  total  work  of  buckling  is  as  given 
in  equation  ( 86) • 

The  total  potential  energy  is 

V  =  L'J  -  AW 
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vhen  small  imperfections  of  both  the  twist  and  the  lateral  deflection  are 
considered;  the  total  deformations  are  the  sum  of  the  initial  and  elastic 
components 


where 

w  =  elastic  lateral  deflection 

e 

w  =  initial  lateral  deflection 

o 

cpg  =  elastic  angle  of  rotation 

^o  =  initial  angle  of  rotation. 

Now  equation  (87  )  becomes 


For  the  simply  supported  case  in  question,  the  total  deformations  may  be 
represented  by 
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and  the  initial  imperfections  may  be  chosen  as 


%  -  1  Pn 


COS 


nTTx 


n=l 

00 


(91) 


-O’l 


_  mTTx 
B  cos  -=- 
m  L 


m=l 


20 


It  has  been  shown  that  the  resulting  critical  load  is  within  1.5  percent 
of  the  exact  solution  when  only  the  first  term  of  equation  (9°)  is  used. 
Therefore,  this  simplification  is  utilized  here  also.  Moreover,  the  same 
simplification  is  applied  to  the  initial  deformations.  These  first  terms 
of  equations  (90)  and  (91)  are  substituted  into  equation  (89).  This  total 
potential  must  be  minimized  with  respect  to  the  unknown  amplitudes  of  the 
total  deflections .  Then 

yield 


and 


where 


P  LB1  '  AlJ  "  1  •  >A  ■  “lj-  K 

(92) 

2EI  l£f  •  LB1  -  Aj  -  Hflj  .  K 

(93) 

K  =  0.351 

This  factor  is  the  result  of  performing  the  prescribed  integrations. 

If  the  amplitudes  of  the  elastic  lateral  deflection  and  of  angle  of  elastic 
rotation  are  denoted  by  6  and  $,  respectively,  and  substitution  for  the  total 
deformation  amplitudes  is  made  in  equations  (92)  and  (93),  the  following 
relationships  are  obtained: 


PLK6=2GJ§ 

2E!  j~]2.  6  -  PLK($  +  P)  =  0 


(94) 

(95) 


Equations  (94)  and  (95)  are  algebraic  relationships  coupled  in  6  and  4.  If 
the  expression  for  6,  as  obtained  from  equation  (94),  is  substituted  into 
equation  (95),  it  yields 
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which  is  a  relationship  of  the  Southwell  type. 

A  search  of  the  literature  did  not  reveal  any  existing  experimental  data 
on  this  subject.  Thus,  a  special  series  of  tests  was  made  to  check  the 
process.  The  experimental  setup  less  the  displacement  transducer  is  depicted 
in  Figures  6l  and  62.  The  deformation  mode  of  the  beam  is  clearly  seen  in 
this  picture.  The  concentrated  load  was  applied  via  a  dowel  and  yoke  system 
which  is  readily  seen  in  the  figure.  The  dimensions  of  the  beam  together 
with  the  deformations  under  load  are  given  in  Table  XIV.  The  resulting 
Southwell  Plot  is  given  in  Figure  63.  The  critical  load  deduced  from  this 
curve  is  21.4  pounds  and  is  in  excellent  agreement  with  the  theoretical  value 
of  22.1  pounds.  Additional  tests  made  with  a  small  spring  providing  some  side 
restraint,  show  that  under  these  conditions  there  is  an  increase  in  buckling 
load  just  as  was  demonstrated  by  Eayashi  and  Kihira  for  the  strut  supported  in 
the  same  manner.  Full  details  of  this  test  series  are  given  in  Relerence  4l. 

The  presentation  so  far  has  been  concerned  mainly  with  the  use  of  the  South- 
well  technique  in  interpreting  data  obtained  from  uniform  x-section  members 
tested  as  individual  items.  It  must  not  be  inferred,  however,  that  the 
process  is  not  applicable  to  structures  fabricated  from  such  elements.  The 
first  application  to  a  member  in  a  realistic  structure  appears  to  have  been 
made  by  Lundquist  .  In  1939  he  analytically  formulated  a  method  for  deter¬ 
mining  the  critical  loads  for  the  displacement  variable,  normally  the  mid- 
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Beam  Dimensions 

Thickness 
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in. 

Depth 

d  =  1.953 

in. 

Length 

L  =-  30.69 
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point  lateral  deflection,  used  for  the  Southwell  representation  could  ^ 
equally  well  be,  chosen  as  the  angle  of  end  rotation.  In  i960,  Gregory^  ^ 
applied  the  technique  to  frame  structures  under  torsion  and  bending.  Hd  was 
able  to  use  the  procedure  tg,  develop  a  design  system  for  lattice  girders 
liable  to  lateral  buckling. The  routine  followed  was  to  load  the  memb'qir 
in  question  to  a  load  value  less  than  that  required  to  precipitate  insta¬ 
bility  and,  from  the  load-displacement  relationship  thus  established,  to 
determine  the  effective  initial  imperfection  and  the  theoretical  critical 
moment.  From  these  two  experimentally  established  parameters,  he  then  pred- 
dicted  the  probable  buckling  moment  for  this  real  structure. 

13 

Ariaratnam  theoretically  established  the  validity  of  the  Southwell  type 
linear  relationship  for  struts  having  a  variable  flexural  rigidity  and 
various  end  conditions  and  for  plane  frameworks  subjected  to  both  flexural 
and  torsional  buckling.  In  1963  Stevens  and  Schmidt^  proposed  a  design 
scheme  for  frameworks.  In  this  plan  they  combined  the  Southwell  formulation 
with  an  iteration  procedure.  Thus,  they  were  able  to  predict  the  critical 
loads  in  the  compression  components  of  the  structure. 

It  is  clear  from  the  discussion  already  given  that  the  key  to  the  correlation 
of  ejq>eriment  and  theory  lies,  in  all  cases,  in  the  treatment  of  a  slightly 
imperfect  system.  From  this  presentation  it  might  be  inferred  that  Ayrton 
and  Perry  were  the  first  to  realize  that  initial  and  elastic  deformations 
of  a  strut  could  be  associated  with  the  Euler  load  and  the  actual  load.  This 
is  not  so.  Admittedly,  Euler  himself  appears  to  have  given  consideration  to 
the  initially  straight  column  only.  However,  an  awareness  of  the  influence 
of  initial  curvature  has  existed  since  the  beginning  of  the  nineteenth 
centuiy.  Young, ^  in  his  course  of  Lectures  on  Natural  Philosophy  and  the 
Mechanical  Arts  published  in  1807,  gave  a  theorem  with  regard  to  this  question. 
This  theorem,  theorem  323  in  his  Magnum  Opus,  states: 

"If  a  beam  is  naturally  of  the  form  which  a  prismatic  beam  would  acquire, 
if  it  were  slightly  by  a  longitudinal  force,  calling  its  depth,  b,  its  length, 
e,  the  circumference  of  a  circle  of  which  the  diameter  is  unity,  c,  the  weight 
of  the  modulus  of  elasticity,  m,  the  natural  deviation  from  rectilinear  form, 
d,  and  a  force  applied  at  the  extremities  of  the  axis,  f,  the  total  deviation 
from  the  rectilinear  form  will  be 
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If  Young’s  equation  is  rewritten  in  present-day  language,  we  obtain  the 
formulation 
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•where 


Vj.  =  the  total  deflection 


w  =  the  initial  deflection 

°  El2 

P_  =  Euler  Load  ~ 

E  L2 


P  =  axial  load  applied 


This  is  identical,  vith  the  Southwell  formulation,  equation  ( ll),  provided 
it  is  realized  that  in  his  time  Young  perpetrated  Jacob  Bernoulli's  error 
in  determining  the  expression  for  the  second  moment  ('‘Moment  of  Inertia") 
for  the  cross-sectional  area.  Bernoulli  had  fixed  this  for  a  rectangular 
section  as 
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instead  of  the  presently  acceptable  relationship  of 
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There  is  little  in  the  literature  between  1807  and  1886  to  indicate  that 
engineers  had  taken  note  of  Young's  contribution,  but  it  seems  clear  that 
many  were  aware  of  the  lack  of  agreement  between  the  Euler  theory  for  a 
perfect  body  and  reality.  Smith, ^5  a  contemporary  of  Ayrton  and  Perry, 
was  certainly  cognizant  of  the  importance  of  initial  curvature  and  eccentric 
loading  prior  to  1886.  He  recognized,  about  this  time  too,  that  the  load 
displacement  curves  for  struts  subjected  to  axial  compression  were  hyper¬ 
bolic.  This  is  evidenced  from  a  footnote  to  Ayrton  and  Perry's  report5  and 
also  from  Reference  4-5. 


That  the  leaders  in  the  field  of  engineering  mechanics  of  this  era  recognized 
the  problem  and  were  acquainted  with  the  mathematical  relationship  of 
equation  (l4)  is  clear  from  T.  Claxton  Fidler's  text.  °  In  this  treatise  of 
1887  he  wrote :  "Let  R  denote  the  resilient  force  of  the  ideal  column,  equal 
to  EI2/L2,  and  let  A  denote  the  initial  deflection,  and  6  the  elastic 
deflection,  then 


P  =  R  * 


A  +  6 


and 


6  =  A 


R 


(101) 


Clax Lon  Fidler's  equation  (99)  is  identical  to  Southwell's  equation  (l4). 
Therefore,  the  deflection  will  now  have  a  certain  assignable  value  depending 
on  the  load  P;  and  if  the  load  is  gradually  increased  the  column  will 
exhibit  an  increasing  deflection,  or,  in  other  words,  it  will  always  be  in 
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a  condition  of  stable  equilibrium;  and  looking  at  formula  (99)  it  is  evident 
that  when  the  load  P  is  increased  to  vithin  a  little  of  the  fixed  quantity, 

R,  the  factor  P/R  -  P  vill  be  a  very  large  multiple,  and  a  very  small 
initial  curvature  (a)  vill  then  be  sufficient  to  produce  a  comparatively 
large  deflection  of  the  column. 

Without  doubt,  there  vere  other  derivations  of  this  basic  equation  which 
ve  have  not  been  able  to  trace,  but  few,  if  any,  could  be  more  general  than 
that  of  Westergaard.l°  In  1922,  this  author  made  a  very  detailed  study  of 
the  buckling  of  elastic  strictures .  Where  most  previous  and  contemporary 
analyses  dealt  with  buckling  from  the  point  of  view  of  the  differential 
equations  of  equilibrium,  Westergaard’s  analysis  was  based  upon  the  use  of 
a  generalized  Lagrangian  analysis.  He  made,  of  course,  the  usual  assumptions 
with  regard  to  linearity  of  stress-strain  relationships  and  considered  only 
deflections  which  were  small.  In  this  manner,  he  was  able  to  show  that 
the  relationship  was  by  no  means  restricted  to.  the  simple  column.  Unfortun¬ 
ately,  like  his  predecessors,  Fidler, Smith,  5  and  Young, ^  he  did  not 
recognize  the  potential  of  his  formula  for  correlating  theory  and  experiment. 
Neither,  of  course,  did  other  engineers. 

From  analysis  and  by  reference  to  tests  made  on  virtually  all  practical 
cases  of  tubes,  braced  structures,  and  flat  plates  subjected  to  in-plane 
loadings,  it  has  been  established  in  this  report  that  the  simple  Southwell® 
formula 


is  absolutely  general  and  is  of  great  importance  in  interpreting  experimental 
data  in  all  these  areas. 

Its  equal  applicability  in  the  area  of  shell  structures  is  asserted.  This 
question  is  discussed  in  another  report. The  large  displacement  analysis 
briefly  mentioned  for  flat  plates  under  compression  and  for  the  circular 
disc  with  a  central  hot  spot  is  also  a  method  of  vide  generality  and  is 
discussed  further  in  a  paper  to  be  published  shortly. ^7 

It  is  to  be  emphasized  that  the  Southwell  method  calls  for  diligence  in 
measurement  of  the  displacements  and  for  care  in  interpretation  of  the 
results.  It  can  be  a  dangerous  practice,  for  example,  to  use  a  least- 
squares-fit  method  in  constructing  the  line. 
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